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We discuss the Bloch-state solutions of the stationary 
Gross-Pitaevskii equation and of the Bogoliubov equations 
for a Bose-Einstein condensate in the presence of a one- 
dimensional optical lattice. The results for the compressibil- 
ity, effective mass and velocity of sound are analysed as a func- 
tion of the lattice depth and of the strength of the two-body 
interaction. The band structure of the spectrum of elemen- 
tary excitations is compared with the one exhibited by the 
stationary solutions ("Bloch bands"). Moreover, the numeri- 
cal calculations are compared with the analytic predictions of 
the tight binding approximation. We also discuss the role of 
quantum fluctuations and show that the condensate exhibits 
3D, 2D or ID features depending on the lattice depth and on 
the number of particles occupying each potential well. We 
finally show how, using a local density approximation, our re- 
sults can be applied to study the behaviour of the gas in the 
presence of harmonic trapping. 



I. INTRODUCTION 

Cold atoms in optical lattices exhibit phenomena typ- 
ical of solid state physics like the formation of energy 
bands, Bloch oscillations, and Josephson effects. Many 
of these phenomena have been already the object of ex- 
perimental and theoretical investigation in Bose-Einstein 
condensates. For deep potential wells further important 
effects take place like the transition from the superfluid 
to the Mott insulator phase [1,2]. 

The purpose of this paper is to study some structural 
properties of interacting Bose-Einstein condensed dilute 
gases at T — in the presence of ID periodic potentials 
generated by laser fields (ID optical lattices). Unless 
the confinement in the radial direction is very tight, the 
transition to the insulator phase in ID optical lattices is 
expected to take place only for extremely deep potential 
wells. There is consequently a large range of potential 
depths where the gas can be described as a fully coher- 
ent system, in the framework of the mean field Gross- 
Pitaevskii approach to the order parameter. 

We will explicitly discuss the change in the behaviour 
of the system as a function of the optical potential depth, 
ranging from the uniform gas (absence of optical lattice) 
to the opposite regime of deep wells separated by high 
barrieres (tight binding limit). Special emphasis will be 
given to the role of two-body interactions, whose effects 



will be adressed by varying the average density n. 

An important feature produced by the periodic poten- 
tial is the occurrence of a typical band structure in the 
energy spectra. In this paper we will discuss different 
manifestations of such a band structure, including: 

- the energy per particle £j(k) of stationary Bloch- 
wave configurations consisting in the motion of the whole 
condensate and carrying a current constant in time and 
uniform in space ("Bloch bands"). This energy is natu- 
rally parametrized as a function of the quasimomentum 
k which, together with the band index j (j = 1, 2, ...), is 
the proper quantum number of these states; 

- the chemical potential 



On 



(1) 



of the same stationary configurations, where n is the aver- 
age density of the sample. The chemical potential plays a 
crucial role in the determination of the equation of state 
and emerges as a natural output of the solution of the 
Gross-Pitaevskii equation (see Eq.(10) below); 

- the spectrum Tiujj(q) of the elementary excitations 
("Bogoliubov bands") carrying quasi- momentum q. The 
elementary excitations are small perturbations of the sys- 
tem and in general can be calculated with respect to each 
stationary configuration of quasimomentum k. In this 
paper we will limit the discussion to the elementary ex- 
citations built on top of the groundstate configuration 
(k = 0). They can be determined by solving the lin- 
earized Gross-Pitaevskii equations (see Eqs.(25,26)). 

The three band spectra £j, /ij and hujj represent dif- 
ferent physical quantities and have the same dependence 
on quasimomentum only in the absence of two body in- 
teractions. Due to the periodicity of the problem the 
quasi-momentum can be restricted to the first Brillouin 
zone. Still it is often convenient to consider all values 
of quasi-momentum to emphasize the periodicity of the 
energy spectra in quasi-momentum space. 

In Sect. II, we discuss the equation of state and the 
Bloch bands, assuming uniform confinement in the ra- 
dial direction. Special emphasis is given to the behaviour 
of the compressibility and the effective mass. The com- 
pressibility k is defined by the thermodynamic relation 
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where \x is the chemical potential relative to the ground- 
state solution of the Gross-Pitaevskii equation (/i = 



1 



Hj=i{k — 0)). For systems interacting with repulsive 
forces, the optical trapping reduces the compressibility 
of the system since the effect of repulsion is enhanced 
by the squeezing of the condensate wavefunction in each 
well. The effective mass is defined through the curvature 
of the lowest (j = 1) energy band 
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Here and in the following, we omit the band index j when 
we refer to the lowest band (j = 1). The current flowing 
along the direction of the optical lattice, is fixed by the 
relation 



(4) 



and in the long wavelength limit k — > 0, one finds I — ► 
nk/m*, where we have used the notation m* = m*(k = 
0). Equivalently, one has 



£ ( fc ) fc^O £ ( fc = °) + 
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For small intensities of the laser field the effective mass 
m* approaches the bare value m . Instead, for large in- 
tensities the effective mass is inversely proportional to the 
tunneling rate through the barrier separating neighbour- 
ing potential wells and is strongly enhanced with respect 
to the bare value. 

Based on the knowledge of the compressibility and of 
the effective mass, one can calculate the sound velocity 
according to the thermodynamic relation 
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This quantity fixes the slope of the lowest Bogoliubov 
band at low quasi momenta. 

In Sect. Ill, we discuss the behaviour of the elemen- 
tary excitations. The excitation spectrum (Bogoliubov 
spectrum) is obtained from the solution of the linearized 
time-dependent Gross-Pitaevskii equation and develops 
energy bands hujj(q) periodic in quasi-momentum space. 
From the same solutions one can calculate the excita- 
tion strengths Zj(p) relative to the density operator and 
hence the dynamic structure factor S(p,iv) [3]. Here, p 
and fiuj are the momentum and the energy transferred 
by a weak external probe. Differently from the Bogoli- 
ubov energies huij, the strengths Zj(p), and S(p,oj), are 
not periodic functions of p, putting in clear evidence the 
difference between momentum and quasi-momentum. 

We will often compare our numerical results with the 
predictions of the so called tight binding approximation, 
which is reached when the intensity of the laser field gen- 
erating the optical lattice is so high that only the overlap 
between the wavefunctions of nearest-neighbour conden- 
sates plays a role. In the tight binding limit, the lowest 



Bloch and Bogoliubov bands take analytic forms that will 
be discussed explicitly. 

In Sect. IV, we study the effect of the lattice on the 
quantum depletion of the condensate and comment on 
the validity of Gross-Pitaevskii and Bogoliubov theory. 
Depending on the parameters of the problem, we distin- 
guish between different configurations which show 3D, 
2D or ID features. 

Finally, in Sect.V we include the presence of an addi- 
tional harmonic trapping. This is achieved by employing 
a local density approximation to treat the inhomogeneity 
of the density profile due to the harmonic confinement. 
Important applications concern the frequencies of collec- 
tive oscillations. 



II. COMPRESSIBILITY AND EFFECTIVE MASS 

The aim of this section is to calculate the compress- 
ibility and the effective mass as a function of the average 
density of the gas and of the depth of the optical poten- 
tial. Using these quantities, we will also calculate the 
velocity of sound. 

We consider the following geometry: along the z direc- 
tion the atoms feel the periodic potential created by two 
counterpropagating laser fields, while the system is uni- 
form in the transverse direction. The important length, 
momentum and energy scales of the problem are: the lat- 
tice spacing d = Tr/k op t related to the wavevector fc opt of 
the laser field in the lattice direction, the Bragg momen- 
tum qs — frfcopt = Tiw/d which identifies the edge of the 
Brillouin zone and the recoil energy E R = h 2 n 2 /2md 2 , 
corresponding to the energy gained by an atom at rest 
by absorbing a lattice photon. 

The trapping potential generated by the optical field 
can be simply written in the form 



V(z) = 8 E R sin 2 (H) 



(7) 



where s is a dimcnsionless parameter which denotes the 
lattice depth in units of the recoil energy Er. 

The parameter characterizing the role of interactions 
in the system is gn, defined as the two-body coupling 
constant g = A7rh 2 a/m times the 3D average density n. 
Here a is the s-wave scattering length which will be al- 
ways assumed to be positive. Typical values of the ratio 
gn/ Er used in recent experiments [4-7] range from 0.02 
to 1. 

Since the density is uniform in the transverse direction, 
the transverse degrees of freedom decouple from the axial 
ones and the stationary Gross-Pitaevskii equation can be 
reduced to a ID equation of the form 
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<p(z) = 
H<p(z) . 



(8) 
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where the order parameter ip is normalized according to 

S%\v{z)\Hz = \. 

In spite of its non-linearity, Eq.(8) permits solutions in 
the form of Bloch waves 



ikz/K ~ 



(9) 



where k is the quasimomentum, j is the band index and 
the function ipjk(z) is periodic with period d. Notice 
that Eq.(9) does not exhaust all the possible station- 
ary solutions of the Gross-Pitaevskii equation [8]. The 
Gross-Pitaevskii equation (8), rewritten in terms of the 
functions ifijk(z), reads 



+gnd\tf jk {z) 



<Pjk{z) = fj,j{k)ipjk{z). (10) 



From the solution of Eq.(10) one gets the functions 
<f>jk{z) and the corresponding chemical potentials fJ,j(k). 

The energy per particle £j(k) can be calculated using 
the expression 
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(p jk {z)dz. (11) 



and differs from the chemical potential Hj(k). In fact by 
multiplying Eq.(10) by Cp* k and integrating, one finds the 
expression 



fd/2 

(k)= < 
J-d/2 



<P*ik( Z ) 
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-s E R sin 2 (z) + gnd\ip jk (z)\ 2 ] ip jk (z)dz. 
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for the chemical potential which coincides with Eq.(ll) 
only in absence of the interaction term. In general, [ij 
and £j are related to each other by Eq.(l). 

The solution of Eq.(10) for k = and j = 1 gives the 
ground state of the system. This state corresponds to a 
condensate at rest in the frame of the optical lattice. In- 
stead the solutions of Eq.(10) with k ^ 0, describe states 
of the system where all the atoms, occupying the same 
single-particle wavefunction, move together with respect 
to the optical potential giving rise to the constant cur- 
rent (4). Experimentally such states can be created by 
turning on adiabatically the intensity of a lattice moving 
at fixed velocity [4]. In this way, it is possible to map 
higher Brillouin zones onto higher bands. 

Results for the Bloch bands (11) are shown in Fig.l 
for s = 5 and gn/E R — 0, 0.1 and 0.5. The effect of in- 
teractions for these parameter values can be hardly dis- 
tinguished in the energy band structure. It can be made 
more evident by plotting the group velocity 



de{k) 
dk 
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as a function of the quasi-momentum. For gn = O.HEr 
one finds a difference in the group velocity of about 30% 
with respect to the non interacting case. The quantity 
plotted in Fig. 1(b) is accessible experimentally through 
Bloch oscillations experiments [9]. In this context, it is 
important to note that the Bloch states (9) with k ^ 
can become energetically or dynamically unstable de- 
pending on the choice of the values for s and gn (See 
for example [10-14]). Note also that for values of gn > s 
one encounters loops ( "swallow tails" ) in the band struc- 
ture at the band edge of the lowest band [12,15] and, at 
even smaller values of gn, at the center of the first ex- 
cited band [12]. For the values of gn we have considered, 
these swallow tails exist only for very small values of the 
lattice depth s. 




FIG. 1. (a) Lowest three Bloch bands for s — 5, gn = 
(solid line), O.IEr (dashed line) and 0.5Er (dash-dotted line). 
The energy of the groundstate (k = 0) has been subtracted; 
(b) Group velocity for the same parameters as in (a) . 
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Let us now focus on the properties of the ground state 
(k = 0, j = 1). By solving numerically Eq.(10) and 
calculating the chemical potential, it is straightforward 
to evaluate the inverse compressibility (2) as a function 
of the relevant parameters of the problem. The results 
are plotted in Fig. 2 as a function of gn for s = 0, 5, 10. 
The case s = is the uniform case, where the equation 
of state is a = gn and = gn. In the presence of the 
optical lattice, we predict a deviation from this linear de- 
pendence on density. One finds an increase of the inverse 
compressibility with s, which is a direct consequence of 
the localization of the wavefunction at the bottom of the 
wells produced by the optical lattice. When the effect of 
two body interactions on the wavefunction is negligible 
one can account for this increase through the simple law 
= g(s)n. The effective coupling constant g(s) de- 
pends only on the lattice depth s and takes the explicit 
form 

fd/2 
J-d/2 

where f gn =o is the groundstate solution of Eq.(8) for 
gn = 0. It is correct to describe k _1 with a linear de- 
pendence on the density only for small interaction pa- 
rameters gn. For higher values of the interactions the 
slope of the curves tends to decrease. This is due to the 
fact that interactions tend to broaden the order parame- 
ter in each well and hence counteract the effect produced 
by the optical lattice. Numerical results for kT x / 1 gn are 
presented in Fig. (3) and compared with the density in- 
dependent quantity g/g (solid line), confirming that in 
general the compressibility does not depend linearly on 
the interaction. However for large s the density depen- 
dence of k^ 1 j gn becomes less and less important and 
the expression K _1 (n, s) = g(s)n becomes applicable in 
a larger range of gn- values. 

2, . . . . 1 



1.5 




0.2 0.4 0.6 0.8 1 

gn/E R 

FIG. 2. Inverse compressibility k^ 1 — ndfi/dn as a func- 
tion of gn/ER for s = (solid line), s — 5 (dashed line) and 
s — 10 (dashed-dotted line) . 

Let us now determine the effective mass by studying 
the low-fc behaviour of the lowest band e(k). The results 



for m* = m*(k = 0) are shown in Fig. 4. For s — > 0, 
the effective mass tends to the bare mass m. Instead for 
large s the effective mass increases strongly due to the de- 
creased tunneling between neighbouring wells of the opti- 
cal potential. The effect of the interactions is to decrease 
the value of m* as a consequence of the broadening of 
the wavefunction caused by the repulsion, which favours 
tunneling, contrasting the effect of the lattice potential. 
In fact, the effective mass is fixed by the tunneling prop- 
erties of the system, which are exponentially sensible to 
the behaviour of the wavefunction in the region of the 
barriers. Then, any small change in the wavefunction 
due to interactions can have a significant effect on the 
effective mass. 
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FIG. 3. n~ 1 /gn for gn = O.IEr (dashed line) and 
gn = 0.5Er (dashed-dotted line) as a function of the lat- 
tice depth s; comparison with the effective coupling constant 
g/g defined in Eq.(14) (solid line). 
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FIG. 4. Effective mass as a function of lattice depth s for 
gn = (solid line), gn = O.IEr (dashed line) and gn — 0.5Er 
(dashed-dotted line). 

The two quantities calculated above, compressibility 
and effective mass, can be used to calculate the sound 
velocity using relation (6) [12,17]. The corresponding re- 
sults are shown in Fig. 5. We find that the sound velocity 
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decreases as the lattice is made deeper. This is due to the 
fact that the increase of the effective mass is more im- 
portant than the decrease of the compressibility k. The 
solid line in Fig. 5 shows that decreasing the interactions 

the sound velocity approaches the law c = ^jjn/ 

where m* n=Q is the effective mass calculated with gn = 0. 
The sound velocity is in principle measurable by studying 
the velocity of a wavepacket propagating in the presence 
of the optical potential. This could be done for example 
by following the experimental procedure used in [18]. Yet 
note that in deep lattices, non linear effects are expected 
to be important also for small amplitude perturbations 
[11]. 




FIG. 5. Sound velocity as a function of the potential depth 
s divided by the sound velocity in the absence of the optical 
potential (s = 0) for gn = 0.02Er (solid line), gn — Q.IEr 
(dashed line) and gn — 0.5Er (dashed-dotted line). 

Most of the results discussed above can be qualita- 
tively understood by working in the so called tight bind- 
ing approximation which becomes more and more accu- 
rate as the intensity of the laser field increases. Within 
the tight binding approximation we can derive analytic 
or semi-analytic expressions for the energy bands and 
consequently for the effective mass, compressibility and 
sound velocity. 

Using Bloch's theorem, we can write the condensate in 
the lowest band as 



Vk {z)=Y,e ikld,h f{z-ld), 



(15) 



where I is the index of the well and / are the Wannier 
functions. The function f(z) is normalised to unity and is 
orthogonal to the functions centered at different sites. In 
general, it depends on two-body interactions and hence 
on the density. Equation (15) holds for any depth of the 
optical potential. However, in the tight binding regime, 
f(z) is a well localized function. This provides important 
simplifications in the calculation of the relevant quanti- 
ties, since only nearest-neighbour overlap integrals have 
to be considered. 



In our derivation of the tight binding results we will 
include also the interaction terms in the Hamiltonian. By 
substituting Eq.(15) into (10) and using definition (11) 
it follows, after some straightforward algebra, that the 
energy per particle takes the simple tight-binding form 



where e 



e(k) 



'kd 
cos I — 
n 



an energy offset, which depends on s and gn but not on 
k, and 5 is the tunneling parameter defined as 



(16) 



f(z)dz is 
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+ V(z) + 2gndf 2 (z) 



f(z-d)dz. (17) 



Using the same approximations, the chemical potential 
takes the form 



H{k) = no 



5ft cos 



(18) 



where /zo 
and 5,i ■ 



.^. + V(z)+gndf(z 



6 



f(z)dz 
To derive 



4gnd J f(z)f(z - d)dz 
Eqs. (16, 18), we have kept terms of the order J f 3 (z)f(z— 
d)dz and neglected terms of the order J f 2 (z)f 2 (z — d)dz, 
which, for localized functions, turn out to be much 
smaller. 

Note that 6 and 8^ depend on density both explic- 
itly and implicitly through the density dependence of / 
(see Eq.(17)). In all situations considered in this work, 
contributions involving df/dn can be safely neglected. 
This approximation allows us to identify the quantity 
5^-5 = -Agnd J f 3 {z)f(x - d) with ndS/dn. 

To check the accuracy of the tight binding approxima- 
tion, one can compare the numerical results for the first 
Bloch band with expression (16). The parameter 5 is re- 
lated to the curvature of the band at k = 0, i.e. to the 
effective mass defined in Eq.(3), through the important 
relation [19] 



(19) 



Hence, we can evaluate (16) by using the numerical re- 
sults for the effective mass discussed above (see Fig. (4)). 
In this way, we automatically include the correct density 
dependence of the tunneling parameter 5. In Fig. 6, we 
compare the first Bloch energy band with its tight bind- 
ing approximation for gn — 0.5-Er and various values 
of s. The comparison shows that, for this value of the 
interactions, the tight binding approximation is already 
quite good for s = 10 and becomes better and better for 
increasing s. 
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FIG. 6. Lowest Bloch band at gn = 0.5Er for different 
values of the potential depths: s = 1 (a), s = 5 (b) and 
s = 10 (c). The solid lines are obtained by evaluating Eq.(ll) 
using the numerical solution of Eq.(10) while the dashed lines 
refer to the tight binding expression (16). The energy of the 
groundstate (k = 0) has been subtracted. 

Using Eq.(18) for k — 0, we can also derive an expres- 
sion for the inverse compressibility, which reads 

kT 1 = gnd J f{z)dz + 8gnd J f(z)f{z - d)dz, (20) 

where contributions due to df ' jdn have been neglected, 
as previously. Since / f 3 (z)f(z — d)dz is normally much 
smaller than f f 4 (z)dz, the on-site contribution to the 
inverse compressibility (first term in (20)) is usually the 
leading term. 

A tight-binding expression for the density-independent 
effective coupling constant g (14) is obtained from (20) by 
replacing / with the Wannier function f gn =o of the non- 
interacting system. Then k" 1 is linear in the density and 
K _1 /n can be identified with g. Neglecting the overlap 
contribution, g takes the form [17] 



9 = gd 



fL=o( z ) dz - 



(21) 



This shows that in the tight binding regime, the effective 
coupling constant g can be safely estimated replacing the 
Bloch state (f gn =o ln Eq.(14) by the Wannier function 

fgn— 0- 

Deep in the tight binding regime, the function f(z) 
can be conveniently approximated by a gaussian / = 
cxp(— z 2 /2cr 2 )/7r 1 / 4 v / a : , where the width a is chosen such 
as to satisfy the equation 



d 3 1 7r 7T 3 , 1 gn 
tt 6 a 6 d d 6 



TTd 2 1 



2 E R V 2 tt 2 a 2 



= (22) 



accounting for the anharmonicities of 0(z 4 ) of the poten- 
tial wells and for the broadening effect of repulsive two 
body interactions. Within the gaussian approximation, 
the inverse compressibility (20) can be rewritten in the 
simplified form 



gnd 
2¥er 



(23) 



where we have neglected the contribution arising from 
the overlap of neighbouring wavefunctions. If we neglect 
the interaction term in Eq. (22), we obtain = gn, 
with g = gd/V2^cr and a = s" 1 / 4 ^ + l/4y/s)d/n. For 
s = 10, gn = 0.5-Er the approximation (23) differs from 
the exact value of k _1 by less than 1%. 

Note that even though the gaussian approximation is 
useful in estimating the compressibility at high s, it can 
not be employed to calculate the effective mass or, equiv- 
alently, the tunneling parameter 5, which requires a more 
accurate description of the tails of the function /. 



III. ELEMENTARY EXCITATIONS 

In this section we study the spectrum of elementary 
excitations. We will focus our attention on the excita- 
tions relative to the ground state, i.e. to the solution of 
Eq.(10) with j = 1 (lowest band) and k = 0. To this 
aim we look for solutions of the time-dependent Gross- 
Pitaevskii equation of the form 



<p(z,t) = e 



— p-i/it/h 



ip(z) + + vUz^M*] , (24) 



where Uj q and Vj q describe a small perturbation with 
respect to the groundstate condensate <p = (fj=i,k=o- 
At first order in the perturbations, the time-dependent 
Gross-Pitaevskii equation yields the Bogoliubov equa- 
tions 



2m 



+ s E R sin 2 {-j^j - M + 2gnd\(p\ 2 



u jq (z) + 



+ sE R sm 2 {-j ) - n + 2gnd\ip 



d 

+gnd(p 2 v jq (z) = TiWj(q)Uj q (z), (25) 

2 f !TZ \ .. , o„„^l..|2 _ 



v jq (z) + 



+gnd(p* Uj q (z) = —tuvj(q)vj q (z). (26) 

The solutions Uj q and Vj q are Bloch waves (uj q = 
exp(iqz /h)iij q (z) where Uj q is periodic with period d and 
analogously for vj q ). They are labeled by their band 
index j and their quasimomcntum q belonging to the 
first Brillouin zone. Hence, also the Bogoliubov spec- 
trum u>j(q) exhibits a band structure [20,21]. 
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FIG. 7. Bogoliubov bands in the first Brillouin zone for 
s = 1, gn = (solid line) and gn — 0.5E R (dash-dotted line). 
Note that for such a small potential, the gap between second 
and third band is still very small. 

Similarities and differences with respect to the well- 
known Bogoliubov spectrum in the uniform case (s — 
0) are immediate. As in the uniform case, interactions 
make the compressibility finite, giving rise to a phononic 
regime for long wavelength excitations (q — > 0) in the 
lowest band. In high bands the spectrum of excitations 
instead resembles the Bloch dispersion (see Eq.(ll)). The 
differences are due to the fact that in the presence of the 
optical lattice the Bogoliubov spectrum develops a band 
structure. As a consequence, the dispersion is periodic 
as a function of quasimomentum and different bands are 
separated by an energy gap. In particular, the phononic 
regime present at q = is repeated at every even multiple 
of the Bragg momentum qs ■ Moreover, the lattice period 
d emerges as an additional physical length scale. 

In Fig. 7 we compare the Bologoliubov bands at s = 1 
for gn = and gn = 0.5Er. In the interacting case, 
one notices the appearance of the phononic regime in 
the lowest band, while higher bands differ from the non- 
interacting ones mainly by an energy shift. 

In Fig. 8 we compare the lowest Bogoliubov and Bloch 
bands. Clearly, the lowest Bloch band is less affected by 
the presence of interactions than the Bogoliubov band. 
Recall that the Bogoliubov band gives the energy of the 
elementary excitations while the Bloch band gives the 
energy per particle of an excitation involving the whole 
condensate. 

The solid lines in Fig. 9 show how the lowest Bogoliubov 
band changes when the lattice depth is increased at fixed 
interaction. At s = 1 (Fig. 9a), apart from the formation 
of the energy gap close to q = qs, the curve still resembles 



the dispersion in the uniform case: both the phononic 
linear regime and the quadratic regime are visible. When 
the potential is made deeper (s = 5, 10; Fig. 9b, c), the 
band becomes flatter. As a consequence, the quadratic 
regime disappears and the slope of the phononic regime 
decreases. This reflects the behaviour of the velocity of 
sound (see Fig. 5). 
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FIG. 8. Lowest Bloch band (dashed line) and lowest Bo- 
goliubov band (dash-dotted line) bands for s = 1 and 
gn = 0.5Er compared with the energy band without interac- 
tions for s = 1 (solid line). The groundstate energy has been 
subtracted in the case of the Bloch band (dashed line) and of 
the energy band without interactions (solid line). 
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FIG. 9. Lowest Bogoliubov band at gn — 0.5Er for differ- 
ent values of the potential depths: s = 1 (a), s = 5 (b) and 
s = 10 (c). The solid lines are obtained from the numeri- 
cal solution of Eqs. (25,26) while the dashed lines refer to the 
tight-binding expression (29). 
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Similarly to the Bloch energy and chemical potential 
spectra, also for the Bogoliubov excitation spectrum we 
can obtain an analytic expression in the tight binding 
limit. We write the excitation amplitudes in the form 



u q (z)^U q J2^ ld/h f(z-ld), 
i 

v q (z) = V q ^e igld/h f(z-ld), 



(27) 
(28) 



where / is the same function as in (15). Recall, that 
/ dzf 2 — 1. Using expression (15) for the order parame- 
ter, we find the result 



hu>„ 



(29) 



86 s 

2 ( 5 + 2n— | sin- 
cm, 



+2kt 



for the excitation frequencies (lowest band), where S and 
are the tunneling and inverse compressibility param- 
eters defined respectively in (17) and (20). In the same 
limit, the Bogoliubov amplitudes are 



jj = £ g + fag 
q 2y/fiw q e q 



V q = 



2^fhL. 



(30) 
(31) 



where oj q is given by Eq.(29) and e q = 25sin 2 (qd/1%) 
captures the quasimomentum dependence of the Bloch 
energy (compare with Eq.(16)). Note that in deriving 
Eqs.(30,31), we have imposed the normalization condi- 

f,onj%\u q {z)\ 2 -\v q {z)\ 2 dz = l. 

In expression (29) the density dependence of the spec- 
trum shows up in three different ways: 

- first of all, the parameter 5 depends on interactions as 
shown explicitly in Fig. 4, where the quantity m* oc 1/5 
(see relation (19)) is plotted; 

- second, k _1 has a more general dependence on the 
density than the one accounted for by the linear law gn, 
as shown in Fig. 3; 

- third, a contribution due to the density derivative of 
6 appears. However this term is always small: for small 
interactions one has ndS/dn <C 5 while for larger inter- 
actions the inverse compressibility k^ 1 dominates both 5 
and ndS/dn. However, as shown in [13,14] this term can 
significantly affect the excitation frequency calculated on 
top of a moving condensate. 

Fig. 9 compares the numerical data with the approxi- 
mate expression (29) evaluated using the parameters n^ 1 
and 5 calculated in the previous section. The tunneling 
parameter 5 is obtained from the data for the effective 
mass m* through Eq.(19). As already found for the low- 
est Bloch band, for this value of gn, the agreement with 
the tight binding expression is already good for s = 10. 



It is possible to identify two regimes, where the Bo- 
goliubov spectrum can be described by further simplified 
expressions: 

(I) for very large potential depth, the spectrum is 
dominated by the second term in the square brackets 
of Eq.(29). In fact, 5 — > while kT 1 becomes larger and 
larger as s increases. Hence, the spectrum takes the form 



hu q w V 5k 1 



sin 



(32) 



both for large and small gn. Of course for large gn, 
the proper density-dependence of 5 and k" 1 has to be 
taken into account in evaluating (32). Note that the Bo- 
goliubov band becomes very flat since 5 decreases expo- 
nentially for large lattice depth s. Yet, its height de- 
creases more slowly than the lowest Bloch band (16) 
whose height decreases linearly in 5. We also point out 
that in the regime where Eq.(32) is valid, the Bogoliubov 
amplitudes Eqs. (30,31) become comparable in magnitude 
for all q in the first Brillouin zone, even if the excitation 
spectrum(32) is not linear. This implies that all excita- 
tions in the lowest band acquire quasi-particle character 
and the role of interactions is strongly enhanced. 

(II) for small enough gn, one can neglect the density 
dependence of 5 and use the approximation = gn, 
where g was defined in (14) and takes the form (21) in 
the tight binding regime. This yields 



'2 5 sin 2 



2 So sin 2 



2gn 



(33) 



which was first obtained in [22] (see also [11,23,24]). 
Eq.(33) has a form similar to the well-known Bogoliubov 
spectrum of uniform gases, the energy 2(5 sin 2 (qd/2h) 
replacing the free particle energy q 2 /2m. 

The spectrum of elementary excitations can be mea- 
sured by exposing the system to a weak perturbation 
transferring momentum p and energy Tilo. The response 
of the system is described by the dynamic structure fac- 
tor S(p, lu) which, in the presence of a periodic potential 
takes the form 



(34) 



where Zj (p) are the density excitation strengths relative 
to the j th band and hujj(p) are the corresponding exci- 
tation energies, defined by the solutions of Eqs. (25, 26). 
Note that p, here assumed to be along the optical lat- 
tice (z axis), is not restricted to the first Brillouin zone, 
being the physical momentum transferred to the system 
by the external probe. In this respect, it is important to 
notice that, while the excitation energies %ujj{p) are pe- 
riodic as a function of p, this is not true for the strengths 
Zj(p). Starting from the solution of Eqs. (25) and (26), 
the excitation strengths Zj (p) can be evaluated using the 
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standard prescription of Bogoliubov theory (see for ex- 
ample [25]) 



zap) = 



d/2 



d/2 



[u* q (z) + v* q (z)} e^/ h V (z)dz 



where q belongs to the first Brillouin zone and is fixed 
by the relation q = p + 2£qB with i integer and 
the Bogoliubov amplitudes are normalized according to 

/-d/2 I u j9( z )| 2 — \ v jq( z )\ 2 d z = 1- The dynamic structure 
factor of a Bogoliubov gas in an optical lattice has been 
recently calculated in [3]. ft is found that for an inter- 
acting system the strength towards the first band devel- 
ops an oscillating behaviour as a function of the momen- 
tum transfer, vanishing at even multiples of the Bragg 
momentum due to the presence of a phononic regime. 
Moreover, in the presence of interactions, the strength 
Z\ towards the first bands vanishes as \^5k in the limit 
of very deep lattices, due to the quasi-particle character 
of the excitations in the whole band. Finally, the sup- 
pression of the static structure factor at small momenta, 
due to phononic correlations, is significantly enhanced by 
the presence of the lattice. 



A. Link with the Josephson formalism 

In this section, we show that the results for the excita- 
tion spectrum in the tight binding limit obtained in the 
previous section (see Eq.(29)) can be recovered using a 
different formalism based on the Josephson equations of 
motion. These can be derived starting from the ansatz 
for the condensate 



nj(t) 



,iSl(t) 



(36) 



where for sufficiently deep optical lattices the wavefunc- 
tion /; is localized at site I and extends only over nearest- 
neighbouring sites, ni is the time-dependent average den- 
sity at site I where the average is taken over the l- 
th well, n is the average equilibrium density and Si is 
the phase of the condensate at site I. At equilibrium 
the functions /; coincide with the Wannier functions 
f(z — Id) introduced before. In general, they depend on 
the density at the corresponding site, as indicated in Eq. 
(36), and hence might themselves implicitly depend on 
time. Furthermore, these functions are chosen such that 
f ff,(z;ni')fi(z;ni)dz = for / ^ V and 1 for I = V . 

When excitations are present, the phases Si at differ- 
ent sites will be different from each other, indicating the 
presence of a current. Using the time-dependent Gross- 
Pitaevskii formalism one can derive the following equa- 
tions of motion for the density and phase variables 



ni = 

V =1+1,1-1 



j?U' 



■y/nini, sin(S; - Si>) , 



(37) 



Si 



i'=i+i,i-i 



(35) where m = J ft 



-L^ + V(z)+gnd\ft\ 

the time-dependent tunneling parameters S 1 ' 1 ' and 8 1 ^' 
arc directly related to the overlap between two neigh- 
bouring wavefunctions 



ftdz, while 



2 J dr. 



fi 



h z d[ 

2m 



+ V ft, + 



(39) 



+gni d ft\ft\ 2 ft> + gnyd ft\ft,\ 2 ft> 



Ignid / ft\ft\ 2 ft<dz. 



(40) 



11' 



Sfj,, where S 



Note that at equilibrium S 1 ' 1 ' = 6 and 5' 
and 5^ have been previously defined in Eq.(17) and after 
Eq.(18). At equilibrium they do not depend on the sites 
I and V since the wavefunctions ft are all the same. 

In order to obtain the excitation frequencies, one has 
to linearize Eqs. (37,38) around equilibrium: in Eq.(37) it 
is enough to take the value of 5 1 ' 1 at equilibrium; instead 
in Eq.(38) one has to expand to first order in the 
density fluctuations An; 



as 1 / as 1 / 

<y M + - 5 ^-An J + - 5 ii-Anj/. 



dni 



dnv 



(41) 



This procedure allows us to recover exactly Eq.(29). In- 
stead, if one sets ^Jninv = ni from the beginning one ob- 
tains result (32) for the excitation spectrum. This proves 
that the first term in the brackets of Eq.(29) has its phys- 
ical origin in the quantum pressure, because it arises from 
the difference in population n ; — ny between neighbouring 
sites. 

Note that in the form (37,38), the equations of motion 
differ from commonly used Josephson equations in that 
the quantities S 1 ' 1 and S 1 ^ 1 can be time-dependent. In the 
usual treatment, these quantities are calculated at equi- 
librium and one approximates 6 = 5^. The resulting sim- 
plified equations are equivalent to the discrete nonlinear 
Schrodinger equation used for example in [26] to inves- 
tigate nonlinear phenomena like solitons and breathers. 
An approach equivalent to the Josephson formalism pre- 
sented in this section, based on the ansatz (36), has been 
developed independently in [13], where the approxima- 
tions involved are discussed in detail. 

For a system confined in the radial direction, contain- 
ing N atoms per site, it is also convenient to introduce 
the Josephson energy Ej = N5 and the charging energy 
Ec = 2d[i/dN = 2k~ 1 /N which play an important role 
in the physics of Josephson oscillations. 
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IV. QUANTUM FLUCTUATIONS AND 
DEPLETION OF THE CONDENSATE 

The presence of the optical potential may introduce 
phase fluctuations which reduce the degree of coherence 
of the sample. This effect is known to yield spectacular 
consequences in 3D optical lattices, giving rise to a tran- 
sition from the superfluid to the Mott insulator phases. 
Also in the presence of a ID optical lattice one can predict 
interesting effects. First of all, the quantum depletion of 
the condensate increases as a consequence of the increase 
of the effective coupling constant (14) and of the effec- 
tive mass. Eventually, if the tunneling rate becomes very 
small, the system reduces to a ID chain of Josephson 
junctions with a modification of the behaviour of long 
range order affecting the phase coherence of the system. 

Let us consider the problem in a 3D box, with the op- 
tical lattice oriented along the z-direction (Note that the 
quantum depletion of the condensate has been calculated 
in [23,24] for different geometries). The quantum num- 
bers of the elementary excitations are the band index j 
and the quasi-momcntum q along the z direction and the 
momenta p x and p y in the transverse directions. The 
quantum depletion of the condensate can be calculated 
using the Bogoliubov result 



AJV t , 
N tnt 



1 



tot 



■2. 



j 1,Px,Py 



d/2 

dz 
-d/2 



■.Jdxjdy \v j<q<Px 



•I'u 



(r)| 2 , (42) 



where JV to t denotes the total number of atoms, AJV to t is 
the number of non-condensed particles and we sum over 
all bands j, over the quasi-momenta q in the first Bril- 
louin zone and the momenta of elementary excitations in 
the transverse directions p x ,p y allowed by the periodic 
boundary conditions. Eq. (42) describes correctly the 
depletion when Bogoliubov theory is applicable. 

In the thermodynamic limit, the depletion can be cal- 
culated replacing the sum with an integral in Eq.(42). In 
the uniform case, the main contribution to the depletion 



is given by quasi-particles with q 



Px 



Pv 



(mcY 



Yet the convergence is very slow and the integral is sat- 
urated by momenta much higher than mc [27], where 
the dispersion exhibits the quadratic p 2 /2m behaviour. 
This implies that in the presence of the lattice it is pos- 
sible to calculate the depiction as in the uniform case, 
provided all the quasimomenta relevant for the calcula- 
tion of the depletion lie within the first Brillouin zone. 
This zone should include a region beyond the phononic 
regime where the dispersion goes like q 2 /2m*. This con- 
dition is satisfied if the inequality m*c <C qs (or cquiva- 
lently kT x <C 5), corresponding to weak interactions and 
relatively low values of s, is fulfilled. Under this condi- 
tion, we can replace q 2 /2m — > q 2 /2m* and g — > g and 
the quantum depletion takes the generalized Bogoliubov 
form 



AJV t( 
N tni 




- — ,/ — (a 3 n) V2 , 
3 7r / V m 



(43) 



where we have defined a through the relation g = 
4irh 2 a/m. In this regime, the depletion will not be quan- 
titatively very different from the corresponding one in 
the absence of the lattice. The situation becomes more 
interesting for larger optical potential depth, where the 
lattice is expected to affect the coherence properties of 
the system. 

In the regime of deep optical lattices, one can neglect 
contributions to the depletion from higher bands, because 
high energy excitations are particle-like. We are then 
allowed to restrict the sum in (42) to j = 1. In the 
tight binding limit, one can easily generalize expressions 
(27,28) for the Bogoliubov amplitudes in the lowest band 
to account for transverse excitations 

e i(pxx+p v v)/h 

( r ) = j U q!Px , Py J2e z " ld/h f(z-ld), (44) 



V <l;Px ,Py ( r ) 



i{jp x x+p y y)/h 



v q ,p x ,p v J2 eiqld/n f( z - ld )' ( 45 ) 



where f(z) is the same function as in (15) and L is the 
transverse size of the system. Neglecting for simplicity 
contributions arising from ndS/dn, Eq.(29) can be gen- 
eralized in a straightforward way to 



Tno q w V £ o (P± , q) (so (P± ,q) + 2n- 



(46) 



where £o(p±,q) = p 2 L /2m + 2Ssin 2 (qd/2h) andp^_ = P x + 
p 2 . For the amplitudes C^ 9 ,p x ,p„ and V q , PxiPy we find the 
result 



_ e ±hw 
u q , Px , Pv , v q , Px , Py - 2V j—- 



satisfying the 

!-d/2 dz I dx I d y [Kp*,pJ 2 



normalization 



(47) 



condition 



2 ]=1- 



\ V 1,Px,Py I 

We replace again the sum (42) by an integral. This 
corresponds to considering the thermodynamic limit in 
all the 3 directions. Inserting Eqs.(45,47), the calculation 
can be performed analytically and gives 



AJVt, 
JV< 



tot 



2 M- 

d \gn 



(48) 



where G(b) = l/2-\/6/7r+6/2-arctan(V6)(l+6)/7r and, 
for simplicity we have used the approximation = gn 
for the compressibility. Since in the tight binding regime 
the ratio 5/gn is usually small and becomes smaller and 
smaller with increasing lattice depth, the result (48) con- 
verges to 



AJVt, 
JV t 



tot 



a 
d 



(49) 



Result (49) coincides with the 2D depletion of a disc of 
axial extension d and scattering length a, where the axial 
confinement is so strong that the motion is frozen along 
the z-direction. In this limit the actual 3D system is 
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described as a series of separated 2D discs. It is inter- 
esting to note that the depletion remains finite even if 
the tunneling parameter S — > 0. The reason is that be- 
fore taking the limit S — > we have taken the continuum 
limit in the radial direction, which forces the system to 
be coherent. In this case the ratio Ej/Ec — Nk5/2 be- 
tween the Josephson energy and the charging energy (see 
end of the last section) is large and hence coherence is 
maintained across the whole sample. A different result 
would be obtained by fixing N and considering the limit 
6 — > 0. We point out that the dependence on the inter- 
action strength is stronger in the 2D case than in the 3D 
case, since the depletion scales like a (2D) rather than 
a 3 / 2 (3D). 

From Eq.(48), one recovers the 3D result (43) in the 
limit 6/gn — > oo. This limit is by the way only of aca- 
demic interest, since in the tight binding limit, where 
(48) was derived, the ratio 6/gn becomes large only 
if interactions are vanishingly small (For example, for 
gn = 0.02Er and s = 10 one finds 6/gn w 1.). 

Under the assumption made in the continuum approx- 
imation, the depletion is always very small and is upper- 
bounded by the quantity a/d. If the continuum approx- 
imation in the radial direction is not applicable and it 
is crucial to take into account the discretization of the 
sum over the quantum numbers p x and p y in Eq.(42), 
the results are different. This is the case if the num- 
ber of particles in each well is sufficiently small, or if the 
longitudinal size of the system, fixed by the number of 
wells N w , is sufficiently large. The limiting case takes 
place when the contribution arising from the term with 
Px = P y = is the dominant one in Eq.(42). In this case 
the system exhibits typical ID features and one finds the 
result 



AiVt, 



tot 



where 



¥ cd 



2ttHN 



(50) 



(51) 



Here, N is the number of particles per well. The de- 
pendence on the interaction strength in the ID case is 
stronger than in the 2D and 3D cases, since c = yjgn/m* 
and hence the depletion scales like a 1 / 2 . This should be 
compared with the a and a 3 / 2 dependence in 2D and 3D 
respectively. 

The transition to the ID character of the fluctuations 
can be identified by the condition 



a 
d 



' In 



4/V„ 



(52) 



which, for gn = 0.2E R , N w = 200 and N = 500, is 
predicted to occur around s = 30 where the left and 
right side of the inequality become equal to ~ 4%. 



Result (50,51) is strictly linked to the coherence the- 
ory of ID systems, where the off-diagonal 1-body density 
exhibits the power law decay 



n 



(1) 



(|r-r'|)-|r-r'|- 



(53) 



at large distances. If the exponent v is much smaller than 
1, the coherence survives at large distances and the appli- 
cation of Bogoliubov theory is justified. For a superfluid, 
the value of v is fixed by the hydrodynamic fluctuations 
of the phase and is given, at T — 0, by the expression 
(51) [28,29]. In terms of the Josephson parameters (see 
section III A) one can also write v — \J ' Ec /8tt 2 Ej. One 
can easily check that, unless N is of the order of unity 
or m* is extremely large, the value of v always remains 
very small. When the exponent of the power law takes 
the value v = 0.14, corresponding to Ej = 1.62Ec, the 
ID system is expected to exhibit the Bradley-Doniach 
phase transition to an insulating phase where the 1-body 
density matrix decays exponentially [30]. Note however 
that before this transition is reached the depletion (50) 
becomes large and hence Bogoliubov theory is no longer 
applicable. 

To give an example, we set gn = 0.2Er, N = 200 and 
N w = 500 describing a setting similar to the experiment 
of [5]. Bogoliubov theory predicts a depletion of w 0.6% 
in the absence of the lattice (s = 0). At a lattice depth of 
s = 10 the evaluation of Eq.(42), using the tight binding 
results (45,47), and keeping the sum discrete yields a de- 
pletion of w 1.7%. On the other hand, Eq.(48), obtained 
by replacing the sum in Eq.(42) by an integral, yields 
a depletion of ss 2%, in reasonable agreement with the 
full result w 1.7%. The 2D formular (49) instead yields 
w 2.9% depletion, revealing that the system is not yet 
fully governed by 2D fluctuations. With the same choice 
of parameters, the power law exponent (51) has the value 
v = 0.001 and the ID depletion (50) is predicted to be 
s=s 0.6%, significantly smaller than the full value 1.7%. 
This reveals that the sum (42) is not exhausted by the 
terms with p x = p y = 0. In conclusion, one finds that 
for this particular setting, the character of fluctuations is 
intermediate between 3D and 2D, and still far from ID. 
In particular, from the above estimates it emerges that in 
order to reach the conditions for observing the Bradley- 
Doniach transition one should work at much larger values 
of s. 



V. APPLICATIONS TO HARMONICALLY 
TRAPPED CONDENSATES 

The results obtained in the previous sections can be 
used to describe harmonically trapped condensates in the 
presence of an optical lattice: If the trapped condensate is 
well described by the TF-approximation in the absence 
of the lattice and if the axial size of the condensate is 
much larger than the interwell separation d, then one 
can generalize the local density approximation (LDA) to 
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describe harmonically trapped condensates in a lattice. 
Basically, the idea is to introduce the average density 



ni(rj_) 



ld+d/2 



n(r_i_,z)dz , 



(54) 



ld-d/2 



where n(r±,z) is the microscopic density and I is the 
index of the lattice sites. Within the LDA, the chemical 
potential is given by 



m , 



M = Vo P t(ni(rj_)) + —{uj z I d + u> ± r ± ) . 



(55) 



where n pt{ni(r±)) is the chemical potential calculated 
at the average density ni(r±) in the presence of the op- 
tical potential and u z , u)± are the axial and transverse 
frequencies of the harmonic trap respectively. Eq.(55) 
fixes the radial density profile iii(r±) at the l-th site once 
the value of m or, equivalently, the number of atoms 



Ni 



rRi 
= 2nd / 
Jo 



r±dr±ni(r±) . 



(56) 



occupying the l-th well is known. In Eq.(56), Ri is the 
radial size of the condensate at the l-th site, fixed by 
the value of r± where the density n/(rj_) vanishes. This 
procedure avoids the full calculation of the microscopic 
density n(r±, z). 

When equilibrium is established across the whole sam- 
ple we have /i; = fx for all I . Making use of this fact and 
employing that Ni — N to t we can find the dependence 
of \x on the total number of particles iV to t- This proce- 
dure also yields the well occupation numbers Ni and the 
number of sites occupied at equilibrium. In the simple 
case in which the chemical potential exhibits a linear de- 
pendence on density u opt = gn + const, one obtains for 
the radial density profile [31] 



ni(r A 



Q V 



in 



-iv ± r ± 



(57) 



where the chemical potential, apart from a constant, 
is given by fi = hu}(l5N tot ag/ahog) 2 ' 5 /2 with u> = 
(ujxLUyUJz) 1 ^ 3 , aho — \Jh/muj. The well occupation num- 
bers and transverse radii are given by 



N l= N (l-l 2 /l 2 m ) 2 , 



2 Ii2 \l/2 



Ri =Ro{l- l 2 /l 2 m ) 



(58) 
(59) 



where l m — y / 2u/muj 2 d 2 fixes the number 2l m + 1 of 
occupied sites, N = 157V/16/ m and i?o = \j2njraio\. 
The increase of fi due to the optical lattice (g > g) implies 
an increase of the radii Ri . This effect has been observed 
in the experiment of [6] . 

The LDA-based approach not only permits to calcu- 
late equilibrium properties, but also dynamic features of 
macroscopic type. To this purpose one can generalize the 
hydrodynamic equations of supcrfluids by taking into ac- 
count the effects of the lattice. Also in this case one can 



use the concept of the average density profile ni(rj_) as 
defined in Eq.(54). Furthermore, one can replace the dis- 
crete index I by the continous variable z = Id. In this 
way, one can define a smoothed macroscopic average den- 
sity profile n,M(r±,z). The dynamics we are interested 
in then consists of the evolution of the macroscopic den- 
sity tim and of the macroscopic superfluid velocity field 
v whose component along the direction of the lattice is 
defined by the average 



I'z = "PT 



1 

D 

HI 

m D 



fD/2 




/ dz 




l-D/2 


m 



(S(D/2) S{-D/2)) . 



(60) 



Here S is the phase of the order parameter and D is a 
length longer than d, but small compared to the size of 
the system as well as to the wavelength of the collective 
oscillations. On this length scale, the effect of Vho can be 
neglected and all the macroscopic variables are constant. 
Hence the phase difference in (60) can be calculated using 
Bloch states with quasimomentum k (see Eq.(9)) and one 
finds the result v z = k/m [32]. Using this result and 
the fact that the energy change per particle due to the 
presence of a small current is k 2 /2m* — m 2 vl/2m* , the 
use of the LDA yields the result 



E= dr 



-n M v x + 



171 m 



2 nMV y + y 

+ e(n M ) 4 



-n M v z 
nuiVho] 



(61) 



for the total energy of the system. Here e(riM) is the 
equilibrium energy per unit volume calculated at the av- 
erage density um in the absence of the harmonic trap 
and both um and v z are now functions of r±, z and t. 

Starting from the functional (61) we can derive the 
hydrodynamic equations 



d_ 

dt 



n M + d x (v x n M ) + d y (v y n M ) + d z ^» z n M j = (62) 



to, m 9 d //to \m ,\ \ „ , . 



+ 



where Vho = m(uj 2 z 2 + uj 2 _r 2 _)/2 is the harmonic po- 
tential and where we have used the relationship fi opt = 
de(nM)/driM- Eqs. (62,63) generalize the hydrodynamic 
equations derived in [17] to situations in which the effec- 
tive mass is density dependent and the chemical poten- 
tial /U op t has a nonlinear density dependence (see section 
(II)). They can be further generalized to account for 
larger condensate velocities [12]. 

The hydrodynamic equations serve to determine the 
frequencies of collective oscillations associated with a 
density dynamics of the form Um{^, t) — njy/(r) + 
e~ luJt 5n(r). Here, 8n(r) denotes a small deviation from 
the equilibrium average density fiu (r) that is associated 
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with a small velocity field v. Hence, it is appropriate to 
linearize Eqs. (62,63) yielding 




0. 



(64) 



This second order equation involves the derivative 
<9^ op t/<9n|- M directly related to the inverse compressibil- 
ity k~ (see Eq.(2)). In general, the solution of Eq.(64) 
should be found numerically because the density depen- 
dence of to* and du opt /dn is not known in an ana- 
lytic form. Yet, for small enough densities the density- 
dependence of m* can be neglected and du opt /dn\- M = g 
as was discussed in Sect. (II). In this case, the frequencies 
of the collective oscillations can be calculated analytically 
[17]. In particular, they can be obtained from the values 
in the absence of the lattice [33] by simply rescaling the 
trapping frequency along z 



(65) 



This result was obtained theoretically for the dipole os- 
cillation in the tight binding regime in [5] and has been 
confirmed experimentally both for the dipole [5] and 
quadrupole oscillation [34]. 

The LDA defined by (55) can also by employed to 
study Bogoliubov excitations occurring on a length scale 
much smaller than the size of the system. Under this con- 
dition, one can define a local Bogoliubov band spectrum 
given by huij(q;nM('r)) where nif(r) is the locally av- 
eraged density introduced above. The Bogoliubov band 
spectrum can be probed by measuring the dynamic struc- 
ture factor S(p,u>), related to the linear response of the 
system to an external perturbation transferring momen- 
tum p and energy Huj. In LDA, the dynamic structure 
factor reads [35] 



(p,u) = Jdrn M {r)S(p,ui;r) 



(66) 



where S(p, u; r) is the dynamic structure factor (34) cal- 
culated at density hm(j) in the absence of harmonic trap- 
ping. As a result of the averaging implied by Eq.(66), the 
dynamic structure factor of the system is not any more 
given by a series of delta-functions as in the absence of 
the harmonic trap, but instead consists of resonances of 
finite width. The validity of the LDA to describe the lin- 
ear response of the trapped condensate in the absence of 
a lattice has been confirmed by Bragg spectroscopy ex- 
periments where the transferred momentum was larger 
than the inverse of the system size and the duration of 
the Bragg pulse was small compared to the inverse of the 
trapping frequencies [36,37]. Experiments and calcula- 
tions which investigate the regime beyond the validity of 
the LDA have been recently carried out [37,38]. 



In the absence of harmonic trapping, the solution of 
the hydrodynamic equations (62) and (63) permits to 
calculate the propagation of sound waves not only in 
condensates at rest, yielding result (6) for the sound ve- 
locity, but also in condensates moving slowly with re- 
spect to the lattice. Let us start from stationary solution 
with macroscopic velocity v z , corresponding to quasimo- 
mentum k = mv z , with k <C qs, and let us consider 
small changes of the velocity field and of the density 
with respect to the stationary values: v z = v z + 5v z , 
riM = nu + Sum- By linearizing Eqs. (62) and (63) and 
looking for solutions oscillating like <~ e —%(qz-ut) ^ onc 
finds at first order in k 



k 

u = c\q\ +q—- 
to* 



(67) 



where 1/to* = d(n/m*)/dn and c is the sound veloc- 
ity in the condensate at rest (see Eq.(6)). The quantity 
to* gives the k — curvature of the lowest chemical 
potential band (see Eq.(12) with j = 1) according to 
a(k) = /i(0) + fc 2 /2m* for k -> 0. Eq.(67) generalizes the 
usual behaviour of the sound velocity in slowly moving 
frames to account for the presence of the optical lattice. 
The significance of to* for the excitation spectrum has 
been pointed out for any optical potential depth for small 
q and any k in [12] and in the tight binding regime for 
any q and any k in [13,14]. 



VI. SUMMARY 

We have studied Bloch-wave solutions of the Gross- 
Pitaevskii equation in the presence of a one-dimensional 
optical lattice. In particular, we have calculated the 
band structure of both stationary ( "Bloch bands" ) and 
time-dependent linearized ("Bogoliubov bands") solu- 
tions. We have discussed these solutions for different 
choices of the lattice depth sEr and the two body in- 
teraction parameter gn. Special attention has been paid 
to the behaviour of the compressibility and of the effec- 
tive mass. We have shown that the compressibility of 
the system is reduced by the presence of the lattice and 
that its inverse is approximately linear in the density for 
low enough gn or high s. In these regimes, the com- 
pressibility and the chemical potential can be expressed 
in terms of an effective coupling constant g > g, which 
accounts for the squeezing of the condensate wavefunc- 
tion in each well. Concerning the effective mass we have 
found that two body interactions give rise to a signifi- 
cant density dependence which decreases its value with 
respect to the prediction for the non-interacting system. 
The compressibility and the effective mass permit to cal- 
culate the sound velocity whose value is found to decrease 
as a function of the lattice depth, reflecting the exponen- 
tial increase of the effective mass. For the tight binding 
regime, we have complemented the numerical results by 
analytic expressions. 
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Concerning the Bogoliubov bands, we have found that 
in a deep lattice the excitations in the lowest band acquire 
a strong quasi-particle character in the whole Brillouin 
zone, characterized by Bogoliubov amplitudes u <~ v. In 
the tight binding regime, analytic expressions for the low- 
est Bogliubov band and the corresponding Bogoliubov 
amplitudes have also been reported. 

In section IV, we have presented results for the quan- 
tum depletion of the condensate. In particular, we have 
found that in a deep lattice the quantum fluctuations of 
the condensate acquire 2D character, reflecting the trans- 
formation of the system into a series of two-dimensional 
discs. As a consequence, the quantum depletion in- 
creases, but remains small, provided the coherence be- 
tween the discs is maintained. If the lattice depth is in- 
creased further, ID quantum fluctuations become impor- 
tant. Estimates of the corresponding effects with realistic 
values of the parameters have been presented. 

Finally, we have demonstrated the use of a local den- 
sity approximation to study macroscopic static and dy- 
namic properties of harmonically trapped systems in the 
presence of an optical lattice. 

In conclusion, Bosc-Einstcin condensates in optical lat- 
tices share important analogies with solid state systems. 
Differences arise due to the presence of two body inter- 
actions giving rise to important new features even in the 
coherent regime where most particles are in the conden- 
sate. The density dependence of the effective mass, as 
well as the distinction between Bloch, chemical poten- 
tial and Bogoliubov bands are some examples discussed 
in this paper. Natural developments of the present work 
concern the study of the dynamics built on top of moving 
condensates where dynamic instabilities can be encoun- 
tered [10-14]. Another important direction is the study 
of nonlinear effects which might sizeably affect the prop- 
agation of sound in the presence of an optical lattice. 
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